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Table 1: Lattice QCD parameters: number of sites, the inverse coupling j3, the clover coefficient c sw , the 
lattice spacing a and the spatial lattice extent L. 



sites 


P 




a [fm] 


L [fm] 


32 3 x 32 


1.83 


1.761 


0.121(2) 


3.87 



1. Introduction 

The //-dibaryon, which has baryon number B = 2 and strangeness S = —2, predicted by 
R. L. Jaffe in 35 years ago [jl|], is one of the most famous candidates of exotic-hadron. This 
prediction was based on two observations: (i) that the quark Pauli exclusion can be completely 
avoided in the flavor-singlet {uuddss) six quark, and (ii) that the one-gluon-exchange interaction 
produce a large attraction for the flavor-singlet six quark [[I], Search for the //-dibaryon is one 
of the most important challenge for theoretical and experimental physics of the strong interaction 
and the quantum chromodynamics (QCD). At present, it is not clear whether there exists the H- 
dibaryon in nature or not. Although deeply bound //-dibaryon with the binding energy Bh > 7 
MeV from the AA threshold has been ruled out by the discovery of the double A hypernuclei, 
^He [||], possibilities of a shallow bound state or a resonance state in this channel still remain [Q]. 

One promising approach to clarify the //-dibaryon is a numerical simulation of QCD on the 
lattice. There was several lattice QCD calculations on the //-dibaryon as reviewed in ref. [^J (see 
also recent works [||, [7[ |L |9]]). However, there was a serious problem in studying dibaryons on 
the lattice: To accommodate two batyons inside the lattice volume, the spatial lattice size L should 
be large enough. Once L becomes large, however, energy levels of two baryons become dense, so 
that isolation of the ground state from the excited states is very difficult (quite a large imaginary- 
time t is required). All the previous works on dibaryons more or less face this issue. While, 
we employed our original breakthrough approach and searched the //-dibaryon, and then found a 



bound //-dibaryon in the flavor SU (3) limit [10, 11]. Because the flavor SU (3) breaking complicate 



calculation, we began with the flavor SU (3) limit. We review our lattice QCD results briefly in the 
next section. In section 3, we estimate effect of the flavor SU (3) breaking on the //-dibaryon. 



2. H-dibaryon from full QCD simulations on the lattice 

As shown in refs. [12, 13, |4|], one can define and extract potential of two-hadron interaction 
from QCD by using the Nambu-Bethe-Salpeter (NBS) wave function measured on the lattice. For 
example, for a single S-wave channel case, the leading order of the derivative expansion of the 
potential can be obtained as 

1 V 2 V(r,t) 4^{r,t) 
^=27^^-W" Ml - M2 (2 " 1} 

where *F(r,f) is a measured NBS wave function, and M\, M2 and /I are measured hadron masses 
and the reduced mass of them respectively. 
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Table 2: Quark hopping parameter K M i s and measured hadron masses: M ps for pseudo-scalar meson octet, 
M vec for vector meson octet, Mt, ar for baryon octet. 





M ps [MeV] 


M vec [MeV] 


M bar [MeV] 




0.13660 


1170.9(7) 


1510.4(0.9) 


2274(2) 


420/4200 


0.13710 


1015.2(6) 


1360.6(1.1) 


2031(2) 


360/3600 


0.13760 


836.5(5) 


1188.9(0.9) 


1749(1) 


480/4800 


0.13800 


672.3(6) 


1027.6(1.0) 


1484(2) 


360/3600 


0.13840 


468.6(7) 


829.2(1.5) 


1161(2) 


720/3600 
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Figure 1: Left: Potential of the flavor-singlet BB channel extracted from QCD at five values of quark mass. 
Right: The ground state of the flavor-singlet BB channel for each quark mass. 



For lattice QCD simulations with dynamical quarks in the flavor SU (3) limit, we have gener- 
ated gauge configurations at five different values of quark masses. We have employed the renor- 
malization group improved Iwasaki gauge action and the non-perturbatively 0(a) improved Wilson 
quark action. Our simulation parameters are summarized in Table [I]. Hadron masses measured on 
each ensemble, together with the quark hopping parameter K uc j s , are given in Table ||. 

The left panel of Fig.[j] shows the flavor-singlet BB potential extracted from our numerical 
simulation of full QCD at the five values of quark mass [11 1. First of all, one sees that the potential 
is entirely attractive and has an attractive core. This feature does not seen in other BB channels and 
hence is characteristic for this flavor singlet BB channel. This lattice QCD result proves that Jaffe's 
and the other quark model prediction to the flavor-singlet BB channel are essentially correct. 

By solving the Schrodinger equation with these potential in the infinite volume, it turned out 
that there is one bound state for each quark mass The right panel of Figjl] shows the energy 
and the root-mean-square (rms) distance of the bound states. These bound states correspond to 
the //-dibaryon predicted by Jaffe. Our lattice QCD calculation shows that a stable //-dibaryon 
certainly exists in the flavor SU (3) limit world with the binding energy of 20-50 MeV for the 
present quark masses. The rms distance y (r 2 ) is a measure of size of the //-dibaryon, which 
can be compared to the point matter rms distance of the deuteron in nature: 3.8 fm. From this 
comparison, one can get feeling of the //-dibaryon: It is much compact than deuteron. 
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Figure 2: Left: The bar-phase-shift of AA and NE. Right: The invariant-mass-spectrum of AA. These are 
calculated with the lattice QCD potentials and the phenomenological baryon masses. 



3. Effect of the flavor SU(3) breaking on the H-dibaryon 

In the real world, the flavor SU (3) symmetry is broken and the //-dibaryon belongs to the 
S = —2, 1 = 0, B = 2 and J p = + sector, instead of the flavor singlet BB channel. There exist three 
coupled BB channels in this sector: AA, AS and EE. In order to study effect of the flavor SU (3) 
breaking on the //-dibaryon, let us take into account the breaking in a phenomenological way. 
In practice, we combine the physical values of baryon masses and the lattice QCD BB potentials 



extracted at a flavor SU (3) limit with the pseudo-scalar meson mass Mps = 469 MeV JTlfJ . This 
calculation is based on the two assumptions: (i) that the major effect of the SU (3) breaking comes 
from the baryon mass splittings, (ii) and that the qualitative features of the hyperon interactions 
remain intact even with the SU (3) breaking. 

By solving the Lippmann-Schwinger equation for T-matrix, 

T afi =v ap + Yyar G (0) T rP^ q(0) = J_ (3 1} 

y E-Hf + ie 

we obtain S-matrix of the spin-singlet S-wave scattering in the coupled channels. The left panel of 
Fig. U shows the obtained bar-phase-shifts 8,- as a function of energy in the center of mass frame. 
We can clearly see a resonance in AA bar-phase-shift at a little below the AS threshold. This 
resonance corresponds to the //-dibaryon. In fact, when we increase a size of the breaking from 
zero to the physical one by hand, we see that the bound state moves to upward and goes above 
the AA threshold for the empirical flavor St/ (3) breaking. If this phenomenological estimation is 
reasonable, the //-dibaryon will be observed in nature as a resonance in AA channel. Similar results 
are reported in refs. [15, 16j 17]. However, to make a definite conclusion on this point from QCD, 



we need to carry out lattice QCD simulations at the physical point and analyses in the coupled 
channels. Study toward this goal is in progress. 

The right panel of Fig. ^| shows the invariant-mass-spectrum of the AA — > AA process cal- 
culated with the above phenomenological treatment of the flavor SU (3) breaking and the S-wave 
dominance assumption. We can see a peak which corresponds to the //-dibaryon. This result 
demonstrates that there is a chance to confirm the existence of the resonant //-dibaryon in nature 
by experiments which count two A's. Actually, there is a report which says that an enhancement 
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of the two A's production has been observed at a little above the AA thresholds in E224 and E522 
experiments at KEK, although statistics is not sufficient for a definite conclusion [Qj. New data 



with high statistics from an upgraded experiment at J-PARC Q18Q as well as the data from heavy 



ion experiments at RHIC and LHC Q19[ ] will shed more lights on the nature of the //-dibaryon in 
near future. 

References 

[1] R. L. Jaffe, Phys. Rev. Lett. 38, 195 (1977) [Erratum-ibid. 38, 617 (1977)]. 

[2] T. Sakai, K. Shimizu and K. Yazaki, Prog. Theor. Phys. Suppl. 137, 121 (2000) 
[arXiv:nucl-th/99 12063]. 

[3] H. Takahashi et ai, Phys. Rev. Lett. 87, 212502 (2001). 

[4] C. J. Yoon et ai, Phys. Rev. C 75, 022201 (2007). 

[5] I. Wetzorke and F. Karsch, Nucl. Phys. Proc. Suppl. 119, 278 (2003) [arXiv:hep-lat/0208029]. 

[6] Z. H. Luo, M. Loan and X. Q. Luo, Mod. Phys. Lett. A 22, 591 (2007) [arXiv:0803.3171 [hep-lat]]. 

[7] S. R. Beane etal. [NPLQCD Collaboration], Phys. Rev. Lett. 106, 162001 (2011) [arXiv:1012.3812 
[hep-lat]]. 

[8] S. R. Beane et al. [NPLQCD Collaboration], Phys. Rev. D 85, 05451 1 (2012) [arXiv:l 109.2889 
[hep-lat]]. 

[9] S. R. Beane, E. Chang, S. D. Cohen, W. Detmold, H. W. Lin, T. C. Luu, K. Orginos and A. Parreno et 
ai, arXiv: 1206.5219 [hep-lat]. 

[10] T. Inoue et al. [HAL QCD Collaboration], Phys. Rev. Lett. 106, 162002 (2011) [arXiv:1012.5928 
[hep-lat]]. 

[11] T. Inoue et ai [HAL QCD Collaboration], Nucl. Phys. A 881, 28 (2012) [arXiv: 1 1 12.5926 [hep-lat]]. 

[12] N. Ishii, S. Aoki and T. Hatsuda, Phys. Rev. Lett. 99, 022001 (2007) [arXiv:nucl-th/061 1096]. 

[13] S. Aoki, T. Hatsuda and N. Ishii, Prog. Theor. Phys. 123, 89 (2010) [arXiv:0909.5585 [hep-lat]]. 

[14] N. Ishii et al. [HAL QCD Collaboration], Phys. Lett. B 712, 437 (2012) [arXiv: 1203.3642 [hep-lat]]. 

[15] S. R. Beane, E. Chang, W. Detmold, B. Joo, H. W. Lin, T. C. Luu, K. Orginos and A. Parreno, Mod. 
Phys. Lett. A 26, 2587 (2011) [arXiv: 1103.2821 [hep-lat]]. 

[16] P. E. Shanahan, A. W. Thomas and R. D. Young, Phys. Rev. Lett. 107, 092004 (2011) 

[arXiv: 1106.2851 [nucl-th]]; A. W. Thomas, P. E. Shanahan and R. D. Young, arXiv:1111.0114 
[nucl-th]. 

[17] J. Haidenbauer and U. -G. Meissner, Phys. Lett. B 706, 100 (201 1) [arXiv: 1 109.3590 [hep-ph]]; 
J. Haidenbauer and U. G. Meissner, Nucl. Phys. A 881, 44 (2012) [arXiv:l 1 1 1.4069 [nucl-th]]. 

[18] J.K.Ahn and K.Imai; J-PARC proposal P42 

[19] N. Shah [for the STAR Collaboration], arXiv:l 1 12.0590 [hep-ex]. 



5 



